Antiferromagnetic ordering induced by paramagnetic depairing in unconventional 

superconductors 
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Antiferromagnetic (AFM) (or spin-density wave) quantum critical fluctuation enhanced just be- 
low Hc2{0) have been often observed in d-wave superconductors with a strong Pauli paramagnetic 
depairing (PD) including CeCoIns. It is shown here that such a tendency of field-induced AFM 
ordering is a consequence of strong PD and should appear particularly in superconductors with a 
gap node along the AFM modulation. Two phenomena seen in CeCoIus, the anomalous vortex 
lattice form factor and the AFM order in the Fulde-Ferrell-Larkin-Ovchinnikov state, are explained 
based on this peculiar PD effect. 

PACS numbers: 



Recently, the heavy-fermion c?-wave paired super- 
conductor CeCoIns with strong paramagnetic depair- 
ing (PD) has been thoroughly studied from the view- 
point of identifying its high field and low temperature 
(HFLT) phase near the zero temperature depairing field 
Hc2{0) with a possible Fulde-Ferrell-Larkin-Ovchinnikov 
(FFLO) state [H, 0]. On the other hand, this material 
also shows an antiferromagnetic (AFM) or equivalently, 
a spin-density wave ordering Q in the HFLT phase and 
transport phenomena suggestive of an AFM quantum 
critical point (QCP) lying near i/c2(0) [4j, |5|. A similar 
AFM fluctuation enhanced near Hc2{0) has also been de- 
tected in other heavy-fermion superconductors 0] and 
cuprates all of which seem to have strong PD and 
a d-wave pairing. A conventional wisdom on this issue 
will be that, in zero field, the nonvanishing supercon- 
ducting (SC) energy gap suppresses AFM ordering and 
thus that the field-induced reduction of the gap leads to 
a recovery of AFM fluctuation. However, it is difficult 
to explain, based on this picture, why an apparent QCP 
is realized not above but only just below Hc2{0) 0, Q in 
those materials. Rather, the fact that the field-induced 
AFM ordering or QCP close to Hc2 (0) is commonly seen 
in superconductors with a d-wave pairing and strong PD 
suggests a common mechanism peculiar to superconduc- 
tivity in finite fields and independent of electronic details 
of those materials such as the band structure. 

In this paper, we point out that, in nodal d-wave su- 
perconductors, a field- induced enhancement of PD tends 
to induce an AFM ordering just below i?c2(0). Although 
relatively weak PD tends to be suppressed by the quasi- 
particle damping effect brought by the AFM fluctuation, 
strong PD rather favor coexistence of a c?-wave super- 
conductivity and an AFM order. Detailed mechanism of 
this field-induced enhancement of AFM fluctuation or or- 
dering below Hc2 depends upon the relative orientation 
between the moment m of the expected AFM order and 
the applied field H : In m || H case, strong PD change 
the sign of the 0(m^|Ap) term in the free energy for any 
pairing symmetry, just like that of its OdAj*) term lead- 
ing to the first order 7Jc2-transition 0, where m = |m| 
and A are order parameters of an expected AFM phase 



and a spin-singlet SC one. In contrast, the field-induced 
AFM ordering in m _L H, which is possibly satisfied in 
CeCoIng in H _L c 3, is a peculiar event to the d-wave 
paired case with the momentum (k) -dependent gap func- 
tion Wk satisfying = — Wk+q and tends to occur irre- 
spective of the presence of the first order iJc2-transition, 
where Q is the wave vector of the commensurate AFM 
modulation. As a consequences of this PD-induced mag- 
netism, two striking phenomena observed in CeCoIns, an 
AFM order Q stabilized by a FFLO spatial modulation 
and an anomalous flux density distribution in the vortex 
lattice [g], wiU be discussed. 

First, we start from a BCS-like electronic Hamiltonian 
for a quasi two-dimensional (2D) material with A 
and m introduced as functionals. By treating A at the 
mean field level, the free energy expressing the two pos- 
sible orderings is written in zero field case as 

J-(H = 0) = ^i|A(r)|2 -rinTr,,,t,™exp[-i7A™/r], 



HA m = ( 51 ^k,"^k'5a,/3Ck,/3 -^[m(q) • S^(q) 

+ h.c.]) +X:(^|m(q)r+A(q)>I't(q)j, (1) 

where \['(q) = - Ek c_k+q/2,tCk+q/2,i, S^^iq) = 
{cri^)a,i3 Z]k4-q.a^k+Q,/j/2, Ck „ crcatcs a quasiparticle 
with spin index a and momentum k, ai, are the Pauli 
matrices, and the positive parameters g and U are the 
attractive and repulsive interaction strengths leading to 
the SC and AFM orderings, respectively. The disper- 
sion Ek, measured from the chemical potential, satisfies 
ek — — Ck-i-Q + TcSi where a small parameter Si 

measuring an incommensurability of the AFM order was 
introduced. In the case with a nonzero H of our interest, 
the Zeeman energy /ibH • (cr)^.^ needs to be added to 
e-kSa, 13- Below, we focus on either the case, m _L H or 
m II H by choosing the spin quantization axis along H. 
The orbital field effect will be included later. 

To examine an interplay between the AFM and SC 
orderings below, let us consider the Gaussian AFM 
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FIG. 1: Diagrams describing (a) x^"' and (b) x'^^'i where the 
cross denotes the particle-hole vertex on the AFM fluctuation, 
and a double solid line in (a) and (b) denotes a normal and 
anomalous Green's function, respectively. 
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fluctuation term Tm in the free energy J- ^ J>„ = 
T^f,lndet[[/-i(5q^q, -Xq,q'(f^)], wherc 



Xq.q'(f^) 



dT(T,5t(q;r)5,(q';0))e*^^^ (2) 



with a fixed and 5*1/ (q;T) denotes 5',y(q) at imagi- 
nary time T. For the moment, we focus on the Pauli 
limit with no orbital field effect and with uniform A 

in which Xq,q'(f^) = [x("Hq,f^) +x(""Hq,f^)]'5q.q', aud 

■^m = -r^qf,lnX(q,r!), where X-i(q,f^) = f/'^ - 

X^^Hq.^^) - X^'^'Hq-^^), and x^"^ and x^™^ are ex- 
pressed by Fig.l (a) and (b), respectively. The A- 
dependent term, x. (q, ^) = X^"'(q, f^)-X^"[iq, ^^)|a=o + 
x(™)(q, O), has been studied in literatures [1^ in = 
case. The corresponding result of Fig.l in a nonzero H 
is given, in the Pauli limit, by 



= T 



p e,(T=±l 

E E - 



4+o,g(ep+Q)4,a(ep) 

i?e+fl,^(P + Q)^e,.(p)' 

(~«;pU;p+Q) |A|2 



,(3) 



where e (il) is a fcrmion's (boson's) Matsubara fre- 



quency, 4,<T(ep) = ie 



FIG. 2: (Color online) Possible Xs = lines (dashed curves) in 
t v.s. ft = H/ Hc2{Q) diagram in H ± c || m following (a) from 
the 0(|Ap) term of Xs(A) for 5i = and (b) from the fuU 
X3(A) in the Pauli limit for 5i = 0.44 (lower dashed (green) 
curve) and 0.63 (higher dashed (blue) one), respectively. With 
no AFM phase in > Hc2{T) (dotted curve), a dashed curve 
is the upper limit of a possible uniform AFM phase. In (a), 
the solid curve is a uniform AFM phase boundary following 
from the dashed curve and an assumed Xq in the normal 
state. The lower panel of (a) is the corresponding XQ^{h) at 
t = 0.0075. However, the uniform AFM order parameter m in 
the uniform A is transformed, in the modulated FFLO state 
with A((^), to the modulated one in{C,): In (b), the upper 
limit (lower solid (blue) curve) of the modulated AFM phase 
for 5i = 0.63 lies close to the second order FFLO transition 
[T^ (higher solid (red)) line. The lower panel of (b) is the 
normalized (JJ-'in*^' (ft) at t = 0.075. 



with a gap node along Q where Wp+qWp < : In 
this case, the 0(|Ap) term of x^"HO,0) remains nega- 
tive and becomes -N{0)\A\'^ /[2{fj.BH)'^] in T ^ limit 
with no PD-induced sign change, where A^(0) is the den- 
sity of states in the normal state. Instead, the cor- 
responding x^**"' (0, 0) and thus, Xs are divergent like 
N{0)[\A\/{hbH)? |ln[Max(t, \Si\)]\ in T limit while 
keeping their positive signs, where t = T/Tc- This di- 
vergence is unaffected by including the orbital depairing. 



a^J.BH + Cp, and D^^^{p) 
AwpP. Main features of x^"-* and That is, in the dj.2_y2-wa.ve case with Q = (tt, tt) [13|, the 



X 



(an) 



are seen in their 0(|A| ) terms. In zero field, 
Xs{A) = Xs(0,0) behaves like T-^ in T ^ limit and 
is negative so that the AFM ordering is suppressed by 
superconductivity flo| . 

To explain results in the case of strong PD, let us first 
focus on m II H case in which a — a. For a near-perfect 
nesting, the 0(|A|2) terms of x'''^-X^"^|a=o and at 
|q| = il = 0, take the same form as the coefficient of the 
quartic (OdAI"*)) term of the SC Ginzburg-Landau (GL) 
free energy and thus, change their sign upon cooling [l^. 
Thus, Xs(A) becomes positive for stronger PD, leading to 
a lower J-„i, i.e., an enhancement of the AFM ordering in 
the SC phase. As well as the corresponding PD-induced 
sign-change of the ©(JAI"*) term which leads to the first 
order i/c2-transition the PD-induced positive Xs is 
also unaffected by inclusion of the orbital depairing. 

In m _L H where a = —a, a different type of PD- 
induced AFM ordering occurs in a d-wave pairing case 



AFM order tends to occur upon cooling even in m _L H. 
In contrast, x^^"HOj 0) is also negative in the da;y-wave 
case satisfying WpWp+q > so that the AFM ordering 
is rather suppressed with increasing H. 

A possible AFM phase boundary in m || c _L H case 
defined as the temperature at which Xq^ = [X{0,0)]~^ 
vanishes up to 0(|Ap) is shown in Fig. 2 (a) together with 
the corresponding Xs = curve, which is the upper limit 
of a PD-induced AFM phase to occur and shifts to higher 
temperature for larger \6i\. Here, the orbital depair- 
ing brought by the gauge field A has been incorporated 
through the quasiparticle Green's function Se,cr(ri , r^) 
in real space with its semiclassical replacement [2|, |l5| 
Gs,<7{ri,r2) = Ge.a{ri - r2) exp[i e{hc)'^ J^^ A ■ dl]. In 

Fig. 2(a), we have used a'^^ (Maki parameter in H _L c) 
= 7.1iiBH^f]^{0)/{2TTTc) = 7.8 and the anisotropy 7 (de- 



fined [14 from Ck) = 4.5, where H^^f^ (0) = jH'^^f' (0) 



r(GL), 



c2,\\ 
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7fi.c/(2e^g) is the iTc2 (O)-value in H _L c defined near Tc 
with the coherence length £,q. The AFM phase is lost 
in H > Hc2 due to the discontinuous _ffc2-transition [2] 
in t < 0.2f5. With a larger Xq^ in the normal state, 
the uniform AFM phase boundary in Fig. 2(a) is pushed 
down to t = to reduce to an AFM-QCP. 

The limitation to the 0(|Ap) term of Xs overestimates 
the AFM ordering. In fact, using the full expression of 
Xs, the AFM region for = is limited to an invisibly 
narrow region in the vicinity of iJc2(0)- Nevertheless, 
the PD-induced AF order for nonzero \6i\ also follows 
from the full Xs(A): By substituting |A| obtained from 
the gap equation into eq.(3), the lines on which the full 
Xs(A) vanishes are obtained, in the Pauli limit with uni- 
form A, as the dashed curves in Fig. 2(b). The Xs > re- 
gion becoming wider with increasing \di\ suggests that, as 
in CeCoIns [3|, the PD-induced AFM order near i?c2(0) 
should be incommensurate. Similar results also follow 
from the direct use of the tight-binding model [iGj . 

Interestingly, the AFM ordered region is expanded fur- 
ther by the presence of the FFLO modulation A(C) = 
■\/2cos((7loC), where C is the component parallel to H of 
the coordinate. To demonstrate this, the gradient expan- 
sion [lH will be applied to the O(m^) term of the mean 
field expression of J^m- Up to the lowest order in the 
gradient, its A(C)-dependent part simply becomes 



-1^(2) 



m(C), 



(4) 

where (Xq ^)'^' = dX ^(k, 0)/9/c^|k=o, and we only have 
to examine its sign which, for uniform A, corresponds 
to that of — Xs- For simplicity, by using the qLO-data 
in Ref.[llj, we find that ST^^'' < below the blue 
solid curve in Fig. 2(b), indicating that the FFLO order 
can coexist with the AFM order in most temperature 
range. Close to the FFLO transition (red solid) curve 
on which qlo = Oj i^iC) is found to have the modula- 
tion ~ sin(gLoC)j so that this AFM order is continuously 
lost as the FFLO nodal planes goes away from the system 
[Tt} . Oppositely, the form of to(C) deep in the FFLO state 
may not be examined properly in terms of the present 
gradient expansion and will be reconsidered elsewhere. 

Next, we examine the anomalous flux density distri- 
bution in the vortex lattice of CeCoIns in H || c Q 
as another phenomenon suggestive of an AFM fluctua- 
tion enhanced just below Hc2- The fiux distribution is 
measured by the form factor which is the Fourier 
component of the longitudinal magnetization Mdr) = 
{Bc{r) — H) I (47r) at the shortest reciprocal lattice vector 
K and implies the slope of the flux density in real 



space. Here, Mc(r) — X]k#o -^^c(K)e'^ '" is given by 



5F 



<5Spd(-K) 



(5) 



B„d=0 



where j(K) and i3pd(K) are Fourier components of the 
orbital supercurrent density and a Zeeman held i?pd im- 
posed to deflne Afc, respectively. It is straightforward 



to, based on the approach in Ref.'^l, obtain Mc in the 
weak-coupling model with no AFM fluctuation by fully 
taking account of the PD and the orbital depairing. In 
our numerical calculations, the OdAj"*) contributions to 
Mc are also incorporated. Nevertheless, it is instructive 
to flrst focus on its 0(|A|2) terms. In i? < 77^^j|^^(0), 
the contribution to Mc from the spin part (last term) of 
eq.(5) is given by 



Mc(K)(P<^) . ^^|AP(K)ImV^(i) (\ + i^ 



(6) 



which is negative, where t/j^^^z) is the flrst derivative of 
the digamma function. Thus, the PD tends to enhance 
the magnetic screening far from the vortex core. This 
is one of origins inducing a fleld-induced increase of 
However, such an enhancement of \F\ calculated 
beyond the 0(|Ap) terms in the weak coupling model 
is much weaker in the d-wave case than in the s-wave 
one [IBl, although it is the c?-wave material CeCoIns 
which has clearly shown such an enhanced \F\ This 
has motivated us to see how the PD-induced AFM 
fluctuation is reflected in For this purpose, let us 
flrst consider the m _L H case and start with the Pauli 
limit again in which Mc is determined by the last term 
of eq.(5). By noting that this term can be obtained 
as the flrst derivative of the free energy density with 
respect to , it is found that the contribution to Mc 
from the Gaussian AFM fluctuation is positive and pro- 
portional to ~ O(|A(r)|2)/^|i7^qA:(q,0)/r3, implying 
a suppressed screening due to the AFM fluctuation, for 
weak PD {fiBH <C T), while it is negative and given by 
-2AiBT|A(r)|2(^Bi/)-^|ln[Max(|<5i|,t)]| EoEq^(q,^^), 
suggesting an enhancement due to the AFM fluctua- 
tion of the screening and thus, of \F\, for strong PD 
{hbH » T), respectively. 

In Fig. 3, we show h v.s. |Fp curves obtained by in- 
corporating the orbital depairing for both cases with and 
with no contributions to eq.(5) from the AFM fluctuation 
through Jvi, where \F\ was normalized by that in the GL 
region near Tc- Based on the result in Fig. 2(a), the famil- 
iar phenomenological form of X{q,il), N{0)X{q,fl) = 



l/[m^it,h) + (|q X z\Co 



(N)^2 



was assumed 



in calculating the J>„-contribution to eq.(5), where 

m^{t,K) = t + 1 — h/liQCP, and ^q^-* = 0.6^0j and v'^ 
is the squared Fermi velocity averaged over the Fermi 
surface. The remarkable peak just below Hc2 of the red 
solid curve is a reflection of the afore- mentioned growth of 
-)j-(an) in m _L H. A similar |i^| -enhancement also occurs 
in m II H (dashed) curves with increasing H, reflecting 
the afore- mentioned sign change of Xs due to large PD. 
To the best of our knowledge, the favorable direction of 
the moment m in CeCoIns in H || c is unclear at present. 
If, in CeCoIns, the m _L H components are dominant in 
the AFM fluctuation even in H || c, the \F{t,h)\-da,ta, 
growing with increasing fleld and on further cooling 
is a reflection of the dj.2_y2-wave pairing ^l5| accompa- 
nied by a strong AFM fluctuation with Q || (tt, tt). With 
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FIG. 3: (Color online) Normalized form factor curves 
in H II c obtained in the weak-coupling BCS model with no 
AFM fluctuation (dotted curves), in the case with the AFM 
contribution with m || H (dashed ones), and in the corre- 
sponding m _L H case (solid ones) a,t t = 0.02 (highest (red) 
curves), 0.04 (middle(blue)), and 0.08 (lowest (green)). Here, 



ftqcp ~ 0.942, a^' = a)^"' /'y = 6.5, and the typical quasi- 
particle damping rate -K'^T/[4:kF^Q^\inN{t, h)y^^] due to 2D 
AFM fluctuation was used. 



the same Q, the corresponding \F(t,h)\ in the s-wave 
case would become much lower than the dotted ones, 
contrary to the trend in the weak-coupling results [l8| . 
The curves in Fig. 3 have been obtained by neglecting the 
narrow FFLO region in CeCoIns in H || c. Effects of the 
FFLO modulation on will be reported elsewhere. 

It has been argued in previous studies P, [3, [3 that 
the anomalous SC properties in CeCoIns in high fields 
are consequences of strong PD. The present results in- 
dicate that the AFM fluctuation and order indicated in 
more recent measurements [S, 9] on this material just be- 
low Hc2{0) also stem from PD, and thus that its HFLT 
phase and the behaviors suggestive of an AFM-QCP near 
Hc2{0) commonly seen in d-wave superconductors with 
strong PD including CeCoIus, should be under- 

stood on the same footing. For instance, the result in 



J20[. However, 
20| is of first or- 



Fig.2 that the field-induced AFM ordering is discontin- 
uously bounded by the first order i7c2-transition natu- 
rally explains why, in spite of the AFM ordering just 
below Hc2iP) [3*1, the transport properties in CeCoIus in 
H > Hc2 have suggested a QCP notably below i?c2(0) 
We have also shown that, in c?3.2_j,2-wave supercon- 
ductors, the AFM order is significantly enhanced by the 
periodic potential due to the FFLO modulation parallel 
to H irrespective of the size of wave number (7lo- In 
contrast, other works have assumed the presence of an 
attractive 7r-triplet pairing channel as the_only factor in- 
ducing the AFM ordering below Hc2 [1^ 
the field-induced AFM transition in Ref. 
der in contrast to the observation [H in CeCoIns. The 
present work has shown that an incommensurate AFM 
ordering just below i/c2(0) develops without assuming 
such an occurrence of other pairing state and due only 
to strong PD effect. In fact, the strange impurity effects, 
arising even from a nonmagnetic doping, on the HFLT 
phase of CeCoIns [2l| cannot be explained unless the 
HFLT phase has a spatial modulation of the SC order 
parameter unrelated to the AFM order [2^ . 

In conclusion, an AFM ordering or fluctuation en- 
hanced close to i/c2(0), often seen in unconventional su- 
perconductors, is a direct consequence of strong param- 
agnetic depairing and also of their c?-wave pairing sym- 
metry with a gap node parallel to the AFM modulation. 
The AFM ordering enhanced due to the FFLO modula- 
tion suggests that the HFLT phase in CeCoIns is a co- 
existent state of the AFM and FFLO orders. Discussing 
the AFM-QCP issues in systems @, Q with a continu- 
ous behavior around Hc2 based on the present theory is 
straightforward and will be performed elsewhere. 
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